Astron. Nachr. / AN XXX, No. XX, 1 -|6](201 1) / DOI please setDOI! 



Alleviation of catastrophic quenching in solar dynamo model 
with nonlocal alpha-effect 

L.L. Kitchatinov^^ * and S.V. Olemskoy^ 

^ Institute for Solar- Terrestrial Physics, P.O. Box 291, Irkutsk 664033, Russia 
^ Pulkovo Astronomical Observatory, St. Petersburg 196140, Russia 

The dates of receipt and acceptance should be inserted later 

Key words Sun: magnetic fields - stars: magnetic fields - magnetohydrodynamics (MHD) - turbulence 

The nonlocal alpha-effect of Babcock-Leighton type is not prone to the catastrophic quenching due to conservation of 
magnetic helicity. This is shown with a dynamo model, which jointly applies the nonlocal alpha-effect, the diamagnetic 
pumping, and dynamical equation for the magnetic alpha-effect. The same model shows catastrophic quenching when 
the alpha-effect is changed to its local formulation. The nonlocal model shows the preferred excitation of magnetic fields 
of dipolar symmetry, which oscillate with a period of about ten years and have a toroidal-to-polar fields ratio of about a 
thousand. 



1 Introduction 

This paper suggest a solution for the problem of so-called 
'catastrophic quenching' of the alpha-effect of the mean- 
field dynamo theory. The standard alpha-effect of helical 
turbulent motions (Parker' 1955; Steenbeck, Krause & Rad- 
ler J,966 ) with its standard (algebraic) quenching by a mag- 
netic field (MoffattHMl Rii digerg Ml Roberts & Soward 
119751 Rudiger & Kitchatinov'1993') is only a part of the to- 
tal alpha-effect known in the mean-field theory. Friesh et 
al. d 19751 1 were probably the first to notice that not only 
helical motions but small-scale helical magnetic fields pro- 
duce the alpha-effect. This magnetic contribution was later 
recognised to lead to the catastrophic quenching of the to- 
tal alpha-effect (Vainshtein & Cattaneo'1992; Gruzinov & 
Diamond [19941 Blackman & Brandenburg 2002) . 

The term 'catastrophic' is used in the sense that the ar- 
gument of the quenching function is RmB^/iJ^q, not just 
i3^/i?^q as in the case of a standard algebraic quenching 
(Rni = ilt/ri is the ratio of turbulent Tyt to microscopic 
77 magnetic diffusivity, and Scq — u[.^yis is the en- 
ergy equipartition value of a magnetic field). The magnetic 
Reynolds number, Rm, in astrophysical fluids is normally 
so large that even a very small mean field B suppresses the 
alpha-effect in the case of catastrophic quenching. 

This type of magnetic quenching is related to the conser- 
vation of magnetic helicity. Its origin is, briefly, as follows. 
The large-scale magnetic fields generated by the alpha-ef- 
fect dynamos are helical. As the magnetic helicity is con- 
served, small-scale magnetic fields attain helicity equal in 
amount and opposite in sign to that of large-scale fields. He- 
lical small-scale fields produce their own magnetic alpha- 
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effect that counteracts the acting alpha-effect of whatever 
origin so that the total alpha-effect vanishes. Detailed dis- 
cussions of catastrophic quenching can be found in litera- 
ture (cf., e.g., Brandenburg & Subramanian l2005l) . 

Catastrophic quenching presents a serious problem for 
the cosmic dynamo theory. The currently leading idea in re- 
solving the problem for the sun is the evacuation of small- 
scale magnetic helicity from the solar interior by helicity 
fluxes (Vishniac & Cho |2001t Guerrero, Chatterjee & Bran- 
denburg 1201 Oi l and then from the solar corona by coronal 
mass ejections (Brandenburg |2009] l . However, coronal ejec- 
tions can evacuate only a minor part of magnetic helicity 
(KUem, Rust & Seehafer|20igi. 

This paper suggests another solution related to nonlo- 
cal formulation of the alpha-effect. If the region where the 
toroidal field is concentrated is spatially separated from the 
region where the poloidal field is produced by the alpha- 
effect, the large-scale fields generated will not be helical, 
and the problem of catastrophic quenching does not arise. 
To show this, we use a numerical model of the ail-dynamo 
with a nonlocal alpha-effect of the Babcock-Leighton type 
(cf., e.g., Dikpati & Charbonneau 11999 1. Another impor- 
tant ingredient of our model is the diamagnetic pumping of 
large-scale fields. The pumping ensures that toroidal fields 
are concentrated near the base of the convection zone away 
from the near-top region where the alpha-effect is active. 
The dynamical equation for the magnetic part of the alpha- 
effect is involved. This normally leads to a strong (catas- 
trophic) suppression of the dynamo. We actually find catas- 
trophic quenching when changing to a local formulation of 
the alpha-effect in our model. The nonlocal model, however, 
does not show any sign of catastrophic quenching. With the 
largest R^ = 10"* we can apply, the results of the runs with 
the magnetic alpha-effect included or neglected are prac- 
tically the same. The nonlocal model reproduces the main 
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features of the solar cycle. There are also some disagree- 
ments with observations showing ways for model improve- 
ment. 

2 The model 

2.1 Dynamical equation for a„ 

Our dynamo model is based on the mean-field induction 
equation 

dB 

— = Vx{£ + uxB-riVxB), (1) 

where rj is the microscopic diffusivity, u is the mean veloc- 
ity, and E = {u' x B') is the mean electromotive force. 
The electromotive force in its local formulation, 

£ ^ {u' X B') ^ {a^ + aJB + (2) 

includes the magnetic alpha-effect together with the st- 
andard kinetic . 

The heuristic dynamical equation for the q;„ 

da,, . _ „ ^r]t !' £ ■ B ^ a„ 

Rn 



^cq 



(3) 



dt ' ' ' "£2 

(Kleeorin & Ruzmaikin |1982t Blackman & Brandenburg 
120021) can be formulated without specification of a partic- 
ular form of the mean electromotive force £; i is the cor- 
relation length. For a uniform mean field, however, Eq. (|2|i 
applies and the steady value of the total a + a„ can 

be estimated to be 

(4) 



a 



1 



This equation describes catastrophic quenching of the 
alpha-effect. We shall use the Eq. (O in our model of so- 
lar dynamo to see that the dynamo is strongly suppressed 
when the local formulation (|2]) for the alpha-effect is ap- 
plied. Note that the magnetic alpha-effect is local by nature 
in contrast to the kinetic alpha-effect that allows a nonlocal 
formulation (Blackman & Brandenburg |2002l Brandenburg 
& Kapvla 2007) . 

With a nonlocal formulation for the kinetic alpha-effect, 
the toroidal field may be small in the region where the alpha- 
effect is active. The a„ is also small in this case and catas- 
trophic quenching does not happen. We shall see that this 
possibility can indeed be realised in a dynamo model with a 
nonlocal (kinetic) alpha-effect. 

2.2 Dynamo equations 

Our model accounts for the diamagnetic pumping of mean 
fields with an effective velocity of 
1, 



(5) 



Diamagnetic pumping was predicted analytically by Zel- 
dovich ( 119571 ) and Radler (1968). The diamagnetic effect 
of inhomogeneous turbulence lacks pictorial explanation 
but its existence has been confirmed by direct numerical 



simulations (Brandenburg et al. 1199 6', Dorch & Nordlund 
i2001i Ziegler & Riidiger 120031) and by laboratory exper- 
iment with turbulent liquid Sodium (Spence et al. 120071 ). 
The concentration of magnetic fields at the base of con- 
vection zone by the diamagnetic pumping can be important 
for a dynamo (Riidiger & Brandenburg 1 1 9951 Guerrero & 
de Gouveia Dal Pino l2008t Kitchatinov & Rudigerl lOOSl ). 

With allowance for the diamagnetic pumping, the mean 
electromotive force reads 



/r^ V X (y^B) + a^B + A, 



(6) 



where 77,^, — rj+r/t is total magnetic diffusivity and A. stands 
for the contribution of the kinetic alpha-effect that can be 
written as 



^=y" a{r,r')B{r')d^ 



(7) 



With a — a-^6{r — r'), we recover the local alpha-effect. 
The kernel function a for the nonlocal formulation will be 
specified later. 

We will consider an axisymmetric dynamo in a spheri- 
cal shell. The magnetic field in this case can be written as 
a superposition of its toroidal part B and a poloidal field 
defined with a toroidal potential A: 



B = e^B + V X 



A 



r sm ( 



(8) 



where standard spherical coordinates are used and is the 
azimuthal unit vector 

Normalized variables are used. Time is measured in un- 
its of Rq/tio; 7/0 is the characteristic value of total diffu- 
sivity. The magnetic field is normalized to the field strength 
Bq for which nonlinear effects become essential, and the a- 
parameter - to its characteristic value ap . The poloidal field 
potential is measured in units of ao-Bo^lj/^yo- From now 
on, the same notations are kept for the normalized variables 
as used before for their not normalized counterparts, except 
for the fractional radius x — r j Rq and normalized diffu- 
sivity rj = rj^ /rjo. The normalized equation for the toroidal 
field reads 



dB 
'dt 



V fdf dA dfdA 
X \dx do do dx 
Tj d / 1 a(sin6'B) 
'^de Vsh^ de 
ld_ ( dj^xB) 
X dx \ dx 



where 



V 



(9) 



(10) 



is the dynamo number The ail-approximation is applied to 
neglect the alpha-effect in the toroidal field equation In 
this equation, / is the normalized rotation frequency. 



u = e^r sm6Qf{x, 



(11) 
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The equation for the poloidal field with nonlocal alpha- 
effect is written as 



dA ri 



d f I dA 



dt ^™ de V sin 61 



d f dA 



+ a;sin6'cos6' J a{x,x')B{x' ,9) dx' + a^B, (12) 

Xi 

where a^, — a;sinf?(a„/ao) is the normalised magnetic 
alpha-parameter. The integration in this equation is only in 
the radius with the upper limit x. This qualitatively reflects 
the fact that the nonlocal alpha-effect in some point x is con- 
tributed by the buoyant magnetic loops rising from deeper 
layers (x' < x) of the convection zone and that buoyant 
velocities are almost vertical. 

The equation system is closed with the dynamical equa- 
tion for . We neglect the helicity flux in Eq. (O and use 
the afi-approximation to rewrite this equation in normal- 
ized units: 

2 



dt 



'2fi 



B 



dA 
'dt 



VrjdiVrjxB) dA 



fj dA a(sin OB) 
x^ sin 0~d9 de 



X dx dx 

where = Vt/vo = ^^Rm^ is the normalized turbulent dif- 
fusivity. The magnetic alpha-effect and related catastrophic 
quenching can be switched off by omitting the last term in 
Eq. (fT2l) . Later on, we will compare the results obtained 
with account for the magnetic alpha effect and without it. 

The boundary conditions assume an interface with a su- 
perconductor on the inner boundary of radius Xi, 

d {^xB) 



dx 



= 0, A = for X ~ Xi, 



and pseudo-vacuum conditions on the top, 

dA 

B = 0, -T— ^0 for a; = 1. 
dx 



(14) 



(15) 



2.3 Model design 

For the differential rotation, we use the approximation by 
Belvedere, Kuzanyan & Sokoloff (120001) for helioseismo- 
logical data 

1 ^ ^ 
fix, — cos (2m {^-0))Y. ^16) 

m— n— 

The coefficients Cmn of this equation are given in Table 1 of 
Belvedere et al. (I2000I I. Figure[l]shows the angular velocity 
isolines. 

The kernel function of the nonlocal alpha-effect in the 
poloidal field equation (fT2] i was prescribed as follows 

4>hix')(f)aix) 



a{x, x') 



1 + B^{x',t 
1 



Mx') = - (1 - erf ((x' - Xb)//ib)) , 
(t'aix) = l-{l + eri{{x-Xa)/ha)), 




Fig. 1 Angular velocity isolines for the differential rota- 
tion used in the dynamo model. 

where erf is the error function and B^ in the denominator of 
the first equation accounts for the usual algebraic quench- 
ing of the alpha-effect. We always use x^, = Xi + 2.5/ib 
and = 1 — 2.5/iq to ensure smoothness of the kernel 
functions in the simulation domain. The ft-b-parameter is 
the thickness of the near-bottom region of toroidal magnetic 
fields producing the alpha-effect. The is the thickness of 
the near-surface region where this alpha-effect is produced. 
The nonlocal alpha-effect with the kernel function ( fTTI l is 
very similar to the Babcock-Leighton mechanism for the 
poloidal field production used in the dynamo models of Dur- 
ney (119951) and Dikpati & Charbonneau ( 19991 . 

To compare the results of the simulations for nonlocal 
and local alpha-effects we will need a local formulation of 
this effect, which can be obtained by applying the kernel 
function 

2f}5{x — x') 



a[x,x')= _^_^^2 . 

The diffusivity profile of our model reads 



■q{x) = R„ 



1 + erf 



(18) 



(19) 



The aim of this paper is to demonstrate that the nonlocal 
alpha-effect is not prone to catastrophic quenching. We will 
do this with the following parameter values: x^ — 0.7, 



0.74, K = 0.01, K = 0.02, K 



0.002, and 



(17) 



Rq/1 =10. The dependence on the model parameters will 
be discussed in a separate publication. The profiles of dif- 
fusivity and kernel functions ( [TtT i for this parameters set 
are shown in Fig. |2] The maximum value of the magnetic 
Reynolds number we were able to apply is R,n — 10^, and 
all the results of this paper were obtained with this value. 

The system of dynamo equations (|9]), ( fT2] i and (fT3] l was 
solved numerically with the grid-point method and explicit 
time-stepping. The diamagnetic pumping and low diffusion 
in the near-bottom region leads to a high concentration of 
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Fig. 2 Profiles of the normalised total diffusivity and the 
kernel functions ([TtI i of the nonlocal alpha-effect. 

the magnetic field near the bottom. To resolve the fine near- 
bottom structure, a nonuniform grid over the radius with the 
grid spacing Ax ^ r\^l'^ was applied. The grid over the 
latitude was uniform. All the results of the next Section do 
not depend on numerical resolution. This was checked by 
repeating the runs with doubled resolution. 

Equatorial symmetry was usually not prescribed. The 
field was evolved in time starting from a mixed-parity initial 
field and the solution relaxed eventually to a certain equa- 
torial symmetry. In order to determine the critical dynamo 
numbers for excitation of the dipolar {B{G) = —B{tt — 9)) 
and quadrupolar {B{9) = B{tt — 9)) dynamo modes, addi- 
tional boundary conditions selecting the field mode of cer- 
tain equatorial symmetry were imposed on the equator. 

3 Results and discussion 

Our computations show that catastrophic quenching is 
present for local formulations of the alpha-effect but it does 
not exist for the model with a nonlocal alpha effect. 

Figure[3]shows the results of computations for the local 
alpha-effect of Eq. ( fTsT l. The upper line in this Figure was 
obtained with the a„ put to zero in the poloidal field equa- 
tion ( fT2] i. Saturation of the dynamo in this case is due to 
the usual algebraic quenching of the alpha-effect. The lower 
line represents the computation with allowance for the mag- 
netic alpha, a^,, governed by the dynamical equation ( fT3] ). 
Magnetic energy considerably decreases in this case indi- 
cating catastrophic quenching of the dynamo. 

In the model with the nonlocal alpha-effect of Eq. ( fTTl ), 
saturation of the dynamo is due to the standard algebraic 
quenching only. The results of Fig. |4] for a^^ neglected or 
included are practically the same. When algebraic quench- 
ing is switched off by omitting B^ in the denominator of 
the first equation of ( fT7] i, but a^^ is kept finite, field growth 
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Fig. 3 Time dependencies of the volume-averaged square 
of the toroidal field (B^) for the runs with algebraic = 
0, upper line) and catastrophic (a„ 7^ 0) quenching of the 
local alpha effect of Eq. ( fTsl l. Time is measured in units 
of Rq/tiq. The computations were performed for dynamo 
number D — 5.3 x 10^ slightly above the critical value of 
Dc = 5.0 X 10". 

0.006 f « 
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Fig. 4 Time dependencies of the mean square of the 
toroidal field in the model with nonlocal alpha-effect. The 
two bottom lines show the results of the runs with — 
(full line) and for the complete model (a^^ ^ 0, dotted line). 
These two lines are hard to distinguish by eye. The upper 
line is for the computation neglecting the algebraic alpha- 
quenching. The field growth does not saturate in this case in 
spite of the finite The dynamo number D = 2.2 x 10'' 
is slightly above the critical value of Z^c = 1-9 x 10". 

does not saturate. Clearly, the does not play any role in 
these simulations. 

The reason for the inefficiency of a^j in the nonlocal 
dynamo model can be seen from Fig. |5] which shows the 
magnetic field patterns for several instances of a magnetic 
cycle. The toroidal field is highly concentrated at the bot- 
tom. It is small in the near-top region where the alpha-effect 
is active. In this case, the alpha-effect does not contribute to 
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Fig. 5 Toroidal field isolines (top row) and poloidal filed 
lines (bottom row) for several instances of a magnetic cycle. 
The time of the run in units of Rq/tiq is shown at the top. 
Full (dotted) lines show positive (negative) levels and clock- 
wise (anticlockwise) circulation. The pictures of the upper 
row are rescaled so that the upper (dashed) boundary shows 
the radius of r = O.74i?0 below which the toroidal fields 
are localized. D = 2.2 x 10''. 

the first term on the right of the equation (O and the gener- 
ated large-scale fields are not helical. The small-scale mag- 
netic helicity balancing the helicity of large-scale fields is 
not produced and the is small. 

It should be noted that the non-locality of the alpha- 
effect alone does not guaranty nonoccurrence of catas- 
trophic quenching. If the toroidal field in the near-top region 
of the alpha-effect has same sign and same order of mag- 
nitude as the near-bottom toroidal field that produces this 
alpha-effect, then catastrophic quenching can still happen 
(Brandenburg & Kapyla 2007). It can be imagined, how- 
ever, that if such a distributed toroidal field has opposite 
signs near the top and near the bottom, then catastrophic 
amplification of the magnetic field will happen. 

What excludes the catastrophic quenching in our model 
is the joint effect of non-locality and diamagnetic confine- 
ment of the toroidal field in the near-bottom region. This 
region has been long recognised as a favorable site for the 
solar dynamo (cf., e.g.. Oilman [1992| |. The reason why the 
magnetic field should be concentrated at this site was not 
clear, however. We suggest that the concentration can be ac- 
counted for by diamagnetic pumping. 

The poloidal field in Fig.|5]is also concentrated at the 
bottom. The following consideration shows that somewhere 
inside the sun the poloidal field must be much stronger than 
on the surface. The magnitude of the surface poloidal field is 
about 1-2 Gauss only. If the toroidal field is estimated by its 
magnitude in sunspots, it is about 1000 times stronger than 
the surface poloidal field. The solar differential rotation of 
about 30% can produce in the 1 1 years of the solar cycle 
a toroidal field that is at most 40 times stronger than the 
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Fig. 6 Butterfly diagram of the depth-integrated toroidal 
field B of Eq. ( l20l i (top panel) and surface radial field 
(bottom) for the model with nonlocal alpha-effect. Time is 
shown in units of R^/rj^. D = 2.2 x lO'*. 

poloidal filed (the strong radial shear in the tachocline does 
not change this estimation because the radial field should be 
as much weaker there compared to the meridional field as 
the radial shear is larger than the latitudinal shear). There- 
fore, a poloidal field of several tens Gauss is necessary for 
the production of kilogauss toroidal fields. The toroidal field 
in our dynamo model is about 1000 times stronger than 
the polar field due to the near-bottom concentration of the 
poloidal field. 

Figure |6] shows butterfly diagrams for the surface radial 
and deep toroidal fields. The toroidal field diagram shows 
the isolines of the quantity 

1 

B = sme J (j)h{x)B{x) dx, (20) 

to which the Babcock-Leighton surface alpha-effect is pro- 
portional. The factor sin 6 in Eq. ( |20l ) accounts for the de- 
pendence of the length of toroidal flux tubes on latitude (it is 
supposed that the probability of sunspot production is pro- 
portional to the length of the tube). 

The period of the magnetic cycles of Figures |4] and 
|6] equals 0.75i?Q/?7o, that for the diffusivity of 770 ss 
lO^'^ crn^/s is close to the 1 1 year period of the solar cycle. 
The model with nonlocal alpha-effect does not suffer from 
the old problem of too short magnetic cycles typical of the 
local dynamo models (Fig.O. 

The radial field diagram of Fig. |6]is similar to observa- 
tional diagrams of Stenflo (1988) and Obridko et al. (12006] l. 

The field pattern of Fig. |6] is antisymmetric about the 
equator The critical dynamo number for excitation of the 
dipolar modes, = 1.9 x 10**, is substantially smaller 
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than the critical number — 2.5 x 10^ for the quadrupo- 
lar modes. Accordingly, the runs started from an initial field 
of mixed parity rapidly relaxed to dipolar parity. So clear 
preference for the dipolar modes is typical of the dynamo 
models with relative small magnetic diffusion near the base 
of the convection zone (Chatterjee, Nandy & Choudhuri 
120041) . 

The only clear disagreement of the model with observa- 
tions is the presence of toroidal fields in Fig.|6]on too high 
latitudes. It may be expected that even a slow meridional 
flow in the near-bottom region of low diffusion can influ- 
ence the latitudinal profile of a toroidal field. Allowance for 
the meridional flow is a perspective for model improvement. 
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